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I.  INTRODUCTION 


This  report  describes  a computer  curve-fitting  program  for  experimental, 
monotonically  increasing  stress-strain  data.  In  its  present  state  of  develop- 
ment the  program  is  useful  for  stress-strain  curves  which  have  a knee  and 
strain  hardening.  For  data  from  zero  to  approximately  the  0.2%  offset  yield 
strength,  a modified  form  of  the  Ramberg-Osgood1  equation  is  used.  Data 
beyond  the  yield  in  the  strain-hardening  region  are  approximated  by  a simple 
power  law.  The  computer  program  also  makes  calculations  of  other  stress-strain 
quantities  as  will  be  subsequently  described. 

Ramberg-Osgood  Equation 

Only  three  constants  are  required  to  characterize  the  data  in  the  Ramberg- 
Osgood  equation,  which  has  the  form 

e = S/E  + K(S/E)m  (1) 

where  e = engineering  strain 
S = engineering  stress 

and  the  three  constants  are 

E = Young's  modulus 
K = constant 
m = exponent. 

The  constant  K is  a function  of  Young's  modulus  and  a quantity  designated 
in  Reference  1 as  the  secant  yield  strength  Sj.  This  is  defined  as  the  stress 
value  of  the  intersection  with  the  stress-strain  curve  of  a line  through  the 
origin  of  slope  0.7E.  This  particular  value  of  secant  modulus  results  in  a 
value  of  Si  which  is  close  to  the  0.2%  offset  yield  strength.  A variation  of 
Equation  1 is  used  in  the  computer  program  of  the  form 

e = S/E  + (3/7)  (S^E)  (S/Si)m.  (2) 

Two  of  the  original  Ramberg-Osgood  parameters  have  been  retained  in  the  current 
work.  For  the  third  parameter.  Si  has  been  used  rather  than  K.  As  shown  in 
Reference  1,  Si  and  K are  related  by: 

K = (3/7) (E/Si)m_1 . (3) 

The  choice  to  use  Si  rather  than  K as  the  third  parameter  was  influenced 
by  the  fact  that  Si  can  be  directly  associated  with  a physical  quantity,  a 
yield  strength,  while  K cannot.  It  also  serves  as  the  upper  limit  of  stress 
for  which  the  Ramberg-Osgood  equation  is  applicable. 

Values  of  the  exponent  m are  related  to  the  shape  of  the  knee  of  the  stress- 
strain  curve.  These  ordinarily  range  from  10  to  100:  the  larger  is  the  value  of 

m,  the  sharper  is  the  knee  of  the  curve. 

1.  RAMBERG,  W.,  and  OSGOOD,  W,  R,  Description  of  Stress-Strain  Curves  by  Three  Parameters.  NACA  TN  902,  July  1943. 
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Strain-Hardening  Equation 

Analyses  of  stress-strain  data  from  a number  of  different  materials  have 
shown  that  Equation  2 will  closely  approximate  experimental  values  up  to  Sj,  but 
beyond  this  point  there  is  generally  a poorer  fit.  However,  a simple  power  law 
of  the  form 

S = Aep  (4) 

shows  excellent  agreement  with  test  data  in  the  strain-hardening  region  beyond  a 
stress  value  of  Sj.  The  two  constants  of  Equation  4 are  designated  as  follows: 

A = strain-hardening  coefficient 
p = strain-hardening  exponent. 

One  other  constant  is  required  to  characterize  the  data,  namely,  the  end  stress 
or  strain  value  of  the  monotonically  increasing  region  of  the  stress-strain  curve 
(which  may  or  may  not  be  the  end-point  of  the  test). 

Plastic  Strain  Components 

It  is  assumed  that  the  proportional  limit  stress  is  the  upper  limit  of 
elastic  behavior  as  well  as  the  limit  of  proportionality  between  stress  and  strain. 
This  assumption  is  consistent  with  the  ASTM  definition2  of  the  proportional 
limit  as  the  "highest  stress  for  which  the  offset  is  not  measurable  with  the 
instrument  used."  On  the  basis  of  this  definition,  strains  beyond  the  propor- 
tional limit  can  be  considered  to  consist  of  two  components:  an  elastic  compo- 

nent ee  and  a plastic  component  ep. 

The  elastic  component  is  defined  as 

ee  = S/E.  (5) 

The  plastic  component  is 

ep  = e - S/E.  (6) 

It  is  possible  to  relate  the  applied  stress  to  the  plastic  component  of  strain 
by  fitting  equations  of  the  form  of  Equation  4 to  the  appropriate  data.  The 
values  of  the  coefficients  and  exponents  of  such  equations  can  serve  as  quanti- 
tative measures  of  plastic  behavior.  Effects  of  various  test  environments  and 
varying  thermomechanical  treatments  on  the  plastic  response  of  materials  will  be 
manifested  by  varying  values  of  the  parameters. 

In  a few  materials  one  equation  of  the  form  of  Equation  4 serves  to  describe 
the  plastic  properties  from  the  proportional  limit  to  the  maximum  value  of 
engineering  stress.  In  most  materials  a better  fit  of  the  data  results  when  two 
equations  are  used:  one  equation  from  the  proportional  limit  to  some  intermediate 

2.  Standard  Test  Methods  of  Metallic  Materials.  Designation  E8-69,  1974,  Annual  Book  of  ASTM  Standards,  Part  7,  Philadelphia, 
American  Society  for  Testing  and  Materials,  1974,  p.  578-598. 
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point  in  the  yield  region,  and  a second  equation  to  apply  from  the  intermediate 
point  to  the  maximum  values.  In  the  computer  program  a two-equation  description 
is  used.  The  equations  have  the  form 

s = Aiepq  (7) 

which  applies  from  the  proportional  limit  to  the  intermediate  point  and 

S = A2epr  (8) 

which  applies  from  the  intermediate  point  to  the  maximum  value. 

The  basic  calculations  to  determine  the  constants  for  Equations  7 and  8 are 
first  made  with  the  stress  Si  as  the  intermediate  value.  In  an  optional -to- the- 
user  calculation  other  intermediate  stresses  are  used. 

Computed  Stress-Strain  Properties 

The  two  terms  on  the  right-hand  side  of  Equation  2 are  an  elastic  strain 
component  and  a plastic  strain  component.  Quantities  such  as  the  proportional 
limit  and  various  yield  strengths  can  be  calculated  from  the  plastic  strain  term. 
In  the  computer  program,  a value  of  the  proportional  limit  is  computed  on  the 
assumption  that  at  this  value  of  stress  the  plastic  component  will  have  a value 
of  one  microinch/ inch.  This  strain  is  below  the  threshold  of  most  extensometers 
and  is  consistent  with  the  ASTM  definition.  A similar  calculation  is  made  to 
determine  the  0.2%  offset  yield  stress  in  which  the  plastic  component  is  0.2% 
strain.  Calculations  of  the  corresponding  strain  values  are  also  made. 

Statistical  Evaluation 

The  average  deviation  of  the  analytical  from  the  observed  stress  values, 
defined  as  the  square  root  of  the  average  of  the  squared  deviations,  is  computed 
separately  for  each  of  three  regions  of  the  stress-strain  curve:  elastic,  knee 

region  (from  the  proportional  limit  to  a stress  of  Si),  and  strain  hardening  (for 
stress  values  in  excess  of  Si).  The  deviation  values  are  a quantitative  measure 
of  how  well  the  analytic  expressions  fit  the  observed  experimental  data.  A 
similar  statistical  evaluation  is  made  of  the  plastic  strain  component  equations. 

Abort  Criteria 

Included  in  the  program  are  criteria  to  determine  whether  sufficient  data 
have  been  supplied  for  a valid  curve-fitting  analysis  to  be  made.  If  the  criteria 
are  not  met,  either  the  calculations  are  modified  or  the  program  is  aborted. 

Appl i cations 

Since  it  is  possible  to  characterize  a set  of  engineering  stress-strain  data 
with  the  five  constants  of  Equations  2 and  4 and  a value  of  Smax  or  e^x,  it  be- 
comes unnecessary  to  keep  a readily  accessible  file  of  stress-strain  curves.  In- 
stead, the  six  constants  need  only  be  stored  for  ready  access  since  these  can  be 
used  to  reproduce  the  experimental  data.  A file  of  these  constants  on  punch 
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cards,  tape,  or  other  permanent  storage  means  can  become  a data  bank  of  stress- 
strain  properties.  Accessing  of  the  data  can  be  made  in  any  number  of  ways  depend- 
ing upon  how  the  data  are  filed. 

Comparisons  between  the  properties  of  individual  specimens  can  be  made  by 
comparing  the  values  of  the  corresponding  constants  both  for  the  engineering 
stress-strain  data  and  the  plastic  strain  components.  Such  comparisons  could 
also  reveal  the  effects  of  environmental  conditions  of  test,  heat  treatment,  etc., 
on  the  properties  of  materials. 

The  use  of  an  analytic  formulation  for  stress-strain  data  is  of  value  in 
structural  calculations  involving  elastic-plastic  response.  Historically,  the 
Ramberg-Osgood  formulation  was  developed  for  aircraft  design  using  sheet  materials 
where  various  moduli  were  required  in  the  relations  to  predict  elastic  and  plastic 
buckling. 


II.  METHOD  OF  ANALYSIS 

Separate  analyses  are  made  to  obtain  the  appropriate  constants  for  the  various 
analytic  expressions. 

1.  Elastic  modulus  and  intercepts  of  the  elastic  data  on  the  stress  and 
strain  axes. 

2.  Ramberg-Osgood  constants  Si  and  m of  Equation  2. 

3.  Strain-hardening  constants  A and  p for  the  engineering  stress-strain 
curve  of  Equation  4. 

4.  Strain-hardening  constants  Ai , q,  A2,  and  r for  the  stress-plastic  strain 
curves  of  Equations  7 and  8. 

A schematic  flow  chart  is  shown  in  Figure  1.  Some  specific  details  of  the 
analyses  which  are  not  given  on  the  flow  chart  are  described  in  this  section. 

Elastic  Analysis 

A linear  least-squares  analysis  is  performed  on  all  the  data  below  an  input 
proportional  limit  to  obtain  a value  of  Young's  modulus.  The  intercept  value, 
which  is  also  calculated,  represents  the  intersection  of  the  elastic  line  and  the 
stress  axis  as  shown  in  Figure  2.  It  is  assumed,  however,  that  the  data  should 
pass  through  zero  and  that  initial  effects  in  the  extensometer  in  the  low  loading 
region  of  the  test  are  more  likely  to  cause  a nonzero  intersection  than  effects 
in  the  loading  system.  On  the  basis  of  these  assumptions,  the  strain  intercept 
value,  shown  as  e0  in  Figure  2,  is  calculated  and  then  algebraically  subtracted 
from  all  the  input  strain  values  before  further  analyses  are  performed.  This 
forces  the  analytic  stress-strain  curve  to  pass  through  zero.  The  new  strain 
values  are  designated  as  adjusted  strains. 

After  the  Ramberg-Osgood  analysis  is  completed  a new  proportional  limit  is 
specified  based  upon  the  Ramberg-Osgood  equation.  The  elastic  and  Ramberg-Osgood 


4 


analyses  are  then  repeated.  This  second  iteration  obviates  the  need  for  the  user 
to  enter  a precise  value  of  the  proportional  limit  as  input  to  the  program.  The 
only  requirement  is  that  the  input  proportional  limit  be  an  elastic  value. 

Ramberg-Osgood  Constants 

In  determining  the  values  of  Si  and  m,  methods  originally  proposed  in 
Reference  1 were  adapted  for  the  computer.  As  previously  indicated,  and  shown  in 
Figure  3,  the  value  Si  lies  at  the  intersection  of  a secant  modulus  of  slope 
equal  to  0.7E  with  the  stress-strain  curve.  After  Si  has  been  determined, 

Ramberg  and  Osgood  suggested  that  a second  secant  modulus  of  slope  equal  to  0.85E 
be  drawn  to  intersect  the  stress-strain  curve.  The  point  of  intersection  would 
have  the  values  S2  for  stress  and  e2  for  strain.  Once  these  values  are  deter- 
mined they  can  be  substituted  into  Equation  2 together  with  the  values  of  E and  Si 
and  the  equation  can  be  solved  for  m resulting  in 

m = log (3Si (Ee2  - S2)/7)/log(S2/S1) . (9) 

The  probability  is  quite  small  that  in  a set  of  experimental  stress-strain 
data  there  will  be  individual  values  whose  secant  moduli  are  exactly  0.7E  and 
0. 85E.  It  is  more  likely  that  the  values  Si  and  S2  will  lie  between  pairs  of 
adjacent  experimental  points.  These  pairs  of  points  are  located  by  a search 
procedure  which  will  also  locate  data  whose  secant  moduli  are  exactly  the  required 
values. 

Let  the  experimental  data  be  designated  as  Si  and  e^  where  i varies  from  1 
to  N,  and  where  N is  the  total  number  of  data  points.  Assume  also  that  the  data 
are  entered  into  the  computer  in  monotonically  increasing  order  of  stress  value. 

In  the  search  procedure  for  Si,  the  quantity  ASi  is  defined  as 

ASi  = Si  - 0. 7Eei . (10) 

This  quantity  will  be  positive  for  all  values  of  Si  which  are  less  than  Si 
and  it  will  be  negative  for  all  values  of  Si  which  are  in  excess  of  Si  (ASi  = 0 
for  a data  value  whose  secant  modulus  is  exactly  0.7E).  Assume  now  that  the  value 
Si  lies  between  the  data  points  Sj  and  Sj  + i as  shown  in  Figure  4.  The  quantity 
Sj , computed  according  to  Equation  10,  will  be  positive  while  the  quantity  Sj  + i 
will  be  negative.  In  the  computer  program,  the  quantities  Sj  and  Sj+i  are  found 
for  successive  pairs  of  data  points  or  in  a systematic  search  procedure,  and  after 
each  computation  the  AS  values  are  compared  for  sign  until  a change  in  sign  is 
discovered.  This  change  in  sign  locates  the  points  Sj  and  Sj+i. 

The  constants  for  the  equation  of  a straight  line  connecting  Sj  and  Sj+i  are 
now  found.  The  parameter  Si,  which  is  taken  as  the  intersection  of  this  straight 
line  with  secant  modulus  line  of  slope  0.7E,  can  now  be  calculated.  The  linear 
approximation,  which  is  assumed  for  the  stress-strain  curve  between  Sj  and  Sj+i, 
is  reasonable  if  the  difference  between  the  two  stress  values  is  not  large.  This 
has  generally  been  the  case  with  a number  of  sets  of  data  which  have  been  analyzed. 

A similar  search  procedure  is  used  for  finding  S2.  Here  ASi  is  defined  as 
ASi  = si  - 0.85Eei.  (11) 
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Figure  1.  Flow  chart  of  program  (schematic).  Letter  identifiers  A through  E refer  to  Section  IV 
of  the  text.  Coded  indicators  N,  N1,  N2,  N3,  and  NP  are  defined  at  the  end  of  the  chart. 
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No 


GLOSSARY  OF  INDICATORS 


N Number  of  data  points  in  a set. 
One  data  point  consists  of  a 
stress  $ a strain  value 

N1  Number  of  data  sets  to  be 
analyzed 

N2  Option  number  for  output 
(When  N2=0  only  standard 
output  is  printed) 

N3  Units  used  for  input  data 
N3=l:  ksi  $ micro-in/in 
N3=2:  lb  § in. /in. 

N3=3:  ksi  § in. /in. 

NP  Number  of  input  data 
points  for  S>Sj 


(Figure  1 Continued) 


Stress 


Stress 


Figure  2.  Stress  and  strain  intercepts  of  least-squares 
straight-line  fit  of  elastic  data  (schematic). 


. Strain 


Figure  3.  Stress-strain  curve  showing  secant  moduli 
used  in  Ram  berg- Osgood  analysis  (schematic). 


Figure  4.  Location  of  parameter  S-j  between  Sj  and  Sj+-| 
(schematic).  The  quantity  Sj  is  positive  and  Sj+-|  is  negative. 
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The  strain  value  e2  is  calculated  from 


e2  = S2/0.85E. 


(12) 


The  parameter  m is  then  calculated  from  Equation  9. 

It  is  now  possible  to  calculate  a more  precise  value  of  proportional  limit. 
The  second  term  on  the  right-hand  side  of  Equation  2 is  the  plastic  component  of 
strain  given  as 


If  ep  is  set  equal  to  10-6  inch/inch,  a value  of  S can  be  found  from  Equation 
13  which  can  be  considered  as  the  proportional  limit.  This  is  done  in  the  program 
and  a second  elastic  analysis  is  made  using  the  calculated  Spj  as  the  criterion 
value.  Then  a second  Ramberg-Osgood  analysis  is  performed  and  it  is  the  result 
of  the  second  analyses  which  are  printed  out  in  the  program. 

Strain-Hardening  Constants 

Since  the  strain-hardening  relation.  Equation  4,  is  linear  in  logarithmic 
coordinates,  a linear  least-squares  analysis  is  performed  on  the  strain-hardening 
data  to  obtain  the  parameter  p.  In  order  to  assure  that  Equation  4 will  connect 
with  Equation  2,  it  is  assumed  that  the  point  (Sj,  ej)  will  lie  on  the  curve  of 
Equation  4.  On  the  basis  of  this  assumption,  the  value  of  A is  determined  from 


Plastic  Component  Analysis 

The  plastic  components  of  strain  beyond  the  proportional  limit  are  found  by 
subtracting  the  elastic  component  from  the  total  strain,  as  in  Equation  6.  The 
constants  for  Equations  7 and  8,  describing  the  stress  versus  plastic  component 
behavior,  can  now  be  determined. 

Referring  to  Figure  5,  it  is  assumed  that  the  intermediate  stress  S^nter 
is  equal  to  Sj.  Linear  least-squares  analyses  of  the  logarithms  of  the  data 
from  Spj  to  Si  are  now  performed  to  obtain  the  constants  Ai  and  q of  Equation  7 
and  for  the  data  values  in  excess  of  Si  for  A2  and  r of  Equation  8.  The  four 
plastic  component  parameters  are  printed  as  part  of  the  standard  output  of  the 
program. 

It  is  possible  that  the  parameters  for  Equations  7 and  8 determined  for  a 
lower  intermediate  stress  value  than  Si  will  give  a better  fit  to  the  data.  As 
an  optional  feature  of  the  program,  which  the  user  can  specify  in  the  input,  a 
series  of  calculations  are  made  in  which  the  intermediate  stress  is  varied.  For 
the  first  calculation,  Sj_nter  is  chosen  such  that  the  first  three  data  points  be- 
yond Spi  are  used  in  the  analysis  for  Equation  7 and  the  remainder  of  the  higher 
value  data  points  for  Equation  8.  On  each  subsequent  calculation,  Sjnter  is  in- 
creased so  that  one  additional  data  point  is  included  for  the  Equation  7 analysis 
and  one  fewer  for  that  of  Equation  8.  The  calculations  are  continued  until  the 
value  of  Sinter  becomes  Si . 


ep  = (3/7)(Si/E)(S/Si)®. 


(13) 


A = Si/eiP. 


(14) 
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Statistical  Evaluation  of  Fit 


When  a curve-fitting  procedure  has  been  used  on  experimental  data,  it  is 
appropriate  that  a quantitative  measure  be  made  of  the  goodness-of-fit  of  the  data 
to  the  analytic  approximations.  This  is  done  in  the  program  by  the  'calculation 
of  the  average  deviations  of  the  computed  from  the  observed  stress  values. 

The  individual  deviation  is  defined  as 

Dj_  = Si  - Sci  (15) 

where  S^  = an  observed  stress  value 

Sci  = computed  stress  for  the  adjusted  strain  e^. 

The  average  deviation  Dave  for  a range  of  n data  values  is  defined  as 

Dave  = (EDi^/ro)^.  (16) 

In  the  program,  analytical  stress  values  are  calculated,  and  the  individual 
and  average  deviations  are  found  for  the  data  in  the  following  categories: 

A.  Engineering  Stress-Strain  Curve 

1)  elastic  data  for  values  below  Spi 

2)  stress  values  from  Spj  to  Sj 

3)  stress  values  from  Sj  to  Smax* 


Stress 


Figure  5.  Two-equation  analytical  description  of  plastic  strain 
components  with  a variable  intermediate  boundary  stress  (schematic). 
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B.  Stress-Plastic  Strain  Component 

1)  from  Spi  to  S2  using  Equation  7 

2)  from  S2  to  S,nax  Equation  8 

3)  in  the  optional  calculation  with  variable  Sin2;er 

a)  from  Spi  to  Sinter  using  Equation  7 

b)  from  Sinter  to  Smax  using  Equation  8. 

The  statistical  calculations  for  deviations  of  Equations  7 and  8 wh,en  Sinter 
is  variable  are  made  in  a subroutine  of  the  program. 

Analytical  Stress  Calculations 

Direct  calculations  of  stress  can  be  made  for  given  values  of  engineering 
strain  for  the  elastic  data  and  by  the  use  of  Equation  4 for  the  plastic  data. 
Stress  values  for  the  plastic  strain  components  can  be  calculated  from  Equations 
7 and  8.  However,  Equation  2,  the  Ramberg-Osgood  equation,  cannot  be  solved 
explicitly  for  stress  when  a given  value  of  strain  is  specified.  An  iterative 
numerical  scheme  was  therefore  employed  in  the  program  to  calculate  the  analytic 
stresses. 

In  developing  the  numerical  procedure.  Equation  2 was  transformed  to  an 
equation  in  the  variable  Z = Sc/S2 

Z + (3/7) (Z)ra  - Ee/Si  = 8 (17) 

where  Sc  is  the  analytical  value  to  be  calculated. 

The  variable  Z is  nondimensional  and  has  an  upper  limit  value  of  Zmax  = 1. 

The  lower  limit  value  can  also  be  specified  as  Z0  = Sp2/S2.  A trial  value  of  Z 
in  Equation  17,  designated  as  Zt  for  a given  value  or  strain  ej,  will  not  generally 
cause  the  right-hand  side  to  vanish.  Rather  there  will  be  a residual: 

Zt  + (3/7) Ztra  + Eei/Sx  = Residual.  (18) 

It  is  then  assumed  that  a solution  is  obtained  when  the  trial  value  of  Z leads 
to  a residual  which  is  less  than  10“5.  In  the  program  a formalized  scheme  for 
incrementing  Z iteratively  with  progressively  smaller  increments  is  used  to  obtain 
a solution.  With  this  scheme  and  the  criterion  value  of  the  residual,  a solution 
value  of  Z accurate  to  the  fifth  significant  figure  can  be  obtained  after  approx- 
imately twenty  iterations.  A subroutine  of  the  program  is  used  for  these 
computations. 

Tangent  Modulus 

Differentiation  of  Equation  2 with  respect  to  strain  and  rearranging  terms 
results  in 

dS/de  = Et  = E/(l  + (3m/7) (S/S2)m_1) . (19) 
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Since  for  stress  values  up  to  the  proportional  limit,  = E,  Equation  19  is 
applied  over  the  range  Spj<S<Si . It  should  be  noted  that  the  analytic  value  of 
stress  is  used  in  calculating  E^  from  Equation  19.  Differentiation  of  Equation 
4 results  in 


Et  = ApeP"1  (20) 

which  applies  for  S>Si. 

An  output  listing  is  given  of  values  of  the  tangent  modulus  for  every  adjusted 
strain  value;  also  included  are  values  of  the  calculated  stress. 

Poisson's  Ratio 

In  computing  the  value  of  Poisson's  ratio  beyond  the  proportional  limit,  the 
following  assumptions  are  made. 

1.  The  transverse  strain  consists  of  two  components:  one  contributed  by  the 

longitudinal  elastic  strain  component,  and  a second  by  the  longitudinal  plastic 
component . 

2.  The  total  transverse  strain  can  be  found  from  a law-of-mixtures  addition 
of  the  two  components. 

3.  Incompressibility  applies  to  the  plastic  components  so  that  Poisson's 
ratio  for  the  plastic  components  is  0.5. 

4.  The  material  is  isotropic. 

These  assumptions  lead  to  the  following  definition  of  the  transverse  strain: 

et  = pe  = peee  + 0.5ep  (21) 

where  e^  = transverse  strain 

p = computed  Poisson's  ratio 
e = total  longitudinal  strain 
pe  = elastic  Poisson's  ratio 
ee  = elastic  strain  component,  Equation  5 
ep  = plastic  strain  component.  Equation  6. 

Substituting  Equations  5 and  6 into  Equation  21  and  solving  for  p result  in: 

P = (pe  - 0.5)S/Ee  + 0.5.  (22) 

The  elastic  value  of  Poisson's  ratio  is  part  of  the  input  data  to  the  program. 
In  the  computation  the  adjusted  value  of  strain  and  the  corresponding  analytic 
value  of  stress  are  used  in  Equation  22.  The  calculated  Poisson's  ratio  values 
appear  in  the  same  listing  as  the  tangent  modulus. 
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III.  SUMMARY  OF  INPUT  AND  OUTPUT  DATA 


Specific  details  of  the  input  data  and  their  required  format  are  given  in 
Section  V of  this  report.  Samples  of  the  output  printout  are  given  in  Appendix  I 

Input 


Three  types  of  input  information  are  required: 

1.  Control  information  for  the  program  including  coded  indicators  .for : 
the  number  of  data  sets;  the  number  of  data  points  in  a set;  and  their  units  of 
measurement;  and  the  optional  output  desired. 

2.  Descriptive  information  on  the  specimen  and  the  test. 

3.  Numerical  test  data  and  an  approximate  value  of  the  proportional  limit. 

Load,  rather  than  stress,  may  be  input  to  the  program;  in  this  case  the 
specimen  area  is  required  in  the  input.  If  the  Poisson's  ratio  calculation  is 
chosen  as  an  option,  the  elastic  value  of  the  ratio  is  required  as  input. 

Output 

There  is  a standard  printout  which  always  appears  when  the  program  is  run. 

In  addition  there  are  three  tabular  outputs  which  are  components  of  the  optional 
output  as  follows: 

Listing  1.  Comparison  of  Computed  and  Observed  Stress 
Listing  2.  Tangent  Modulus  and  Poisson's  Ratio 

Listing  3.  Plastic  Component  Parameters  with  a Variable  Intermediate  Stress 

The  user  may  select  to  have  none  of  the  three  listings  printed  or  one  of  the 
following  three  options: 

Option  1.  Listing  1 only 
Option  2.  Listings  1 and  2 
Option  3.  All  three  listings 

Standard  Printout 

The  following  information  is  given: 

1.  Specimen  and  test  descriptive  information. 

2.  Young's  modulus  and  the  intercepts  of  the  elastic  data  on  the  stress  and 
strain  axes. 

3.  Ramberg-Osgood  parameters  Sj  and  m,  proportional  limit  stress  and  strain 
and  0.2%  offset  yield  stress  and  strain. 

4.  Strain-hardening  parameters  A and  p for  the  engineering  stress-strain 
curve. 
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5.  Strain-hardening  parameters  Aj , A2,  q,  and  r for  the  stress-plastic 
strain  component  equations  with  S^nter  = Sj. 

6.  Statistical  average  deviations  of  the  calculated  stresses  for  each  of  the 
analytic  expressions. 

Tabular  Listing  1:  Comparison  of  Computed  and  Observed  Stress 

Individual  values  are  given  of  each  of  the  following  quantities. 

1.  Input  stress,  input  strain,  adjusted  strain,  and  computed  stress  from 
Equation  2 or  Equation  4 

2.  Deviation  of  each  computed  stress  value.  Equation  15 

3.  Plastic  strain  component  values  and  computed  stress  from  Equation  7 or 
Equation  8 with  intermediate  stress  of  S^ 

4.  Deviation  of  each  computed  stress  value 

Tabular  Listing  2:  Tangent  Modulus  and  Poisson's  Ratio 

This  listing  includes  individual  values  of: 

1.  Computed  stress  and  adjusted  strain 

2.  Tangent  modulus  and  ratio  of  tangent  modulus  to  elastic  modulus 

3.  Poisson's  ratio 

Tabular  Listing  3:  Plastic  Component  Parameters  with  Variable  Intermediate  Stress 

Included  are  values  of  the  coefficients  and  exponents  for  the  two-equation 
description  of  plastic  strain  components  for  a range  of  variable  intermediate 
stresses  as  follows: 

1.  Aj  and  q for  Equation  7 

2.  Average  deviation.  Equation  16,  for  stresses  predicted  by  Equation  7 

3.  Intermediate  stress  and  associated  plastic  component  of  strain 

4.  A 2 and  r for  Equation  8 

5.  Average  deviation.  Equation  16,  for  stresses  predicted  by  Equation  8 

IV.  PROGRAM  ABORT  CRITERIA 

Calculations  in  the  program  will  be  completely  or  partially  aborted  when  the 
input  data  do  not  meet  certain  criteria.  In  each  case  when  such  a premature 
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termination  occurs,  a specific  diagnostic  statement  will  appear  in  the  data  print- 
out. The  criteria  and  the  resulting  modification  of  the  computing  procedure,  when 
the  criteria  are  not  met,  are  discussed  below  as  they  are  applied  in  the  program. 
The  tests  of  the  data  are  indicated  by  letter  identifiers  on  the  flow  chart  (see 
Figure  1),  which  correlate  with  the  letter  identifier  of  the  paragraph  titles. 

A.  Order  of  Input  Stress  Data  - The  input  test  data  must  be  arranged  in 
increasing  order  of  stress  magnitude.  If  this  criterion  is  not  met,  the  calcula- 
tions are  terminated. 

B.  Elastic  Data  - At  least  three  data  values  below  the  proportional  limit 
are  required  before  the  least-squares  analysis  for  modulus  and  intercepts  will  be 
performed.  If  this  criterion  is  not  met  on  the  first  or  second  iteration,  the 
calculations  are  terminated. 

C.  Input  Data  to  the  Yield  Stress  - The  stress  value  of  the  last  input  data 
point  must  equal  or  exceed  the  secant  yield  strength  Sj.  If  this  criterion  is 
not  met,  the  Ramberg-Osgood  analysis  and  subsequent  analyses  are  terminated  and 
only  the  elastic  data  are  printed  out. 

D.  Knee-Region  Data  - At  least  two  data  points  are  required  between  and 

S2  as  shown  in  Figure  6.  When  fewer  points  are  present  the  analyses  are  terminated 
as  in  C and  only  the  elastic  data  are  printed  out. 

E.  Strain-Hardening  Data  - There  must  be  at  least  three  data  points  whose 
stress  values  are  equal  to  or  greater  than  Sj . If  this  criterion  is  not  met,  no 
strain-hardening  analysis  of  the  engineering  stress-strain  data  is  performed. 

The  calculations  which  are  performed  and  the  output  which  is  printed  then  depend 
upon  what  optional  data  have  been  requested  from  the  program.  The  data  which 
appear  as  output  when  the  various  printout  options  are  requested  are  as  follows: 


Figure  6.  Stress-strain  curve  with  data  points  in  the 
vicinity  of  Si  and  S2  (schematic). 
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1.  No  optional  tabular  data.  Standard  printout  with  the  following  omitted: 
a)  strain-hardening  data  for  Equation  2 and  b)  plastic  component  parameters  for 
Equation  8. 

2.  Option  1.  Standard  printout  as  above  and  tabular  values  for  all  input 
data  less  than  Si  (E i on  Figure  1). 

3.  Option  2.  Same  printout  as  in  Option  1 plus  tabular  values  of  tangent 
modulus  and  Poisson's  ratio  for  all  input  data  less  than  Si  (E2  on  Figure  1). 

4.  Option  3.  Same  as  Option  2,  since  there  are  insufficient  data  beyond 
Si  to  evaluate  the  parameters  for  Equation  8 (E 3 on  Figure  1). 


V.  USER  INPUT  DATA 

The  program  requires  that  the  data  be  entered  into  the  computer  on  punch 
cards.  Specific  details  of  the  input  formats,  columnar  locations,  etc.,  are 
given  in  a supplement*  to  this  report.  A summary  of  the  information  required 
on  the  data  cards  is  given  below  and  listed  in  Figure  7. 

A typical  data  deck  contains  a deck  header  card  followed  by  sets  of  stress- 
strain  data  cards,  each  set  of  which  consists  of  a set  header  card  followed  by 
individual  data  cards,  each  containing  one  value  of  stress  (or  load)  and  one  value 
of  strain. 

First  Card:  The  deck  header  card  requires  the  number  of  data  sets  to  be 
analyzed  (1  to  99  sets)  entered  in  the  first  three  columns  and  the  date  the  program 
is  run  entered  in  columns  4 through  15. 

Data  Cards:  Each  set  of  data  cards  consists  of  a specimen  header  card 

followed  by  the  individual  data  cards.  The  specimen  header  card  requires  the 
following  information: 


Figure  7.  Data  deck  (schematic). 

Maximum  number  of  data  sets  is  99. 
Maximum  number  of  data  cards  per  set  is  99. 


•PAPIRNO,  R.  Supplement  to  AMMRC  TR  76-12.  Program  Listing,  Glossary  and  Specific  User  Instructions,  available  from  Mechanics 
Research  Laboratory  on  request. 
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Specimen  identifier: 

Material  identifier: 

Test  temperature: 

Test  strain  rate: 

Number  of  data  cards: 

Input  format  coded  indicator: 
Proportional  limit  stress: 
Elastic  Poisson's  ratio: 
Specimen  area: 

Date  of  test: 

Output  option  coded  indicator: 


6 columns 
12  columns 
6 columns 

6 columns 
3 columns 
3 columns 

7 columns 
5 columns 

7 columns  (an  entry  is  only  required  if  the 
specimen  load  is  used  on  the  data  cards) 

12  columns 
3 columns 


The  specimen  data  may  be  entered  for  any  one  set  in  one  of  three  forms.  In 
each  case  the  stress  or  load  value  is  entered  in  the  first  seven  columns  and  the 
strain  in  the  next  seven  columns.  The  three  permissible  combinations  are  as 
follows : 

Stress  in  kilopounds  per  square  inch,  strain  in  microinches  per  inch. 

Load  in  pounds,  strain  in  inches  per  inch. 

Stress  in  kilopounds  per  square  inch,  strain  in  inches  per  inch. 

Specific  information  given  in  the  supplement  is  sufficient  to  allow  use  of  the 
program  by  those  who  have  only  a rudimentary  knowledge  of  computer  programming. 
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APPENDIX  I.  TYPICAL  OUTPUT  WITH  OPTION  3 


Sample  of  the  printout  for  a specimen  of  ESR  4340  steel  tested 
in  compression  at  room  temperature.  The  standard  output  and 
the  three  tables  of  Option  3 are  shown. 


ANALYTIC  APPROXIMATION  OF  STRESS-STRAIN  PROPERTIES 
RAHBERG-OSGOOO  ANO  STRAIN  HARDENING  PARAMETERS 
OBTAINED  FROM  ANALYSES  OF  EXPERIMENT AL  DATA 
PROGRAM  NUMBER  ARPY5  MAY  75 

THE  CALCULATIONS  WERE  PERFORMED  ON  MAY  27.1975 


RAHBERC-CSGOOD  EQUATION  EPSrSIC/E*< 3/7  )C Sl/E I ( SIG/S1 1 •• M 

STRAIN  HARDENINC  LAW  SIC” A ( EPS  I • *P 

SPEC.  NO.  ESRD11  4340  STEEL  TEMFr  2C  DEG.  C STRAIN  RATE-  1.C-C5  PER  SEC  TEST  DATE  DEC.  1974 

ELASTIC  PARAMETERS  ELASTIC  MOOULUS.  E-  29.60  MILLION  PSI 

INTERCEPT  ON  STRESS  AXIS-  11C.  PSI  INTERCEPT  ON  STRAIN  AXIS=  -4.  MICRO-IN/IN 

PARAMETERS  OETERMlNEO  AFTER  ALCEBRAIC  SUBTRACTION  OF  STRAIN  OFFSET  FROM  RAW  STRAIN  DATA 


RAHBERG-OSGOOO  PARAMETERS 

Sir 

185.41 

KSI 

Mr  57.08 

PROPORTIONAL  LIMIT  11  MICRO-IN/IN 

OFFSET)  SPLr 

161.46 

KSI 

EPLr  .005454 

IN/IN 

0.2  PERCENT  YIELD 

SYr 

184.45 

KSI 

EYr  .008231 

IN/IN 

STRAIN-HARDENING  PARAMETERS 

COEFFICIENT  Ar 

.2341*06 

EXPONENT  Pr 

.049433 

AVERAGE  STATISTICAL  DEVIATIONS  OF  COMPUTEO  STRESS  VALUES  FROM  INPUT  STRESS  OATA 
FROM  ZERO  TO  SPL  212.  PSI  •••••  FROM  SPL  TO  SI  276.  PSI  •••••  FROM  SI  TO  SMAX  639.  PS  I 

• •••••••• 


STRAIN-HARDENING  PARAMETERS  FOR  PLASTIC  COMPONENTS  OF  STRAIN 
FOR  EQUATIONS  OF  THE  FORM 

SIGrAl«EPS-PLASTIC)**0  FROM  SPL  TO  SI  ••••  SIG rAZ CE PS -PL AST IC 1 • «R  FROM  SI  TO  SMAX 

Al—  • 2066+06  Qr  .01810  A2=  .2147*06  Rr  .02475 

AVERAGE  STATISTICAL  DEVIATION  OF  CCMPUTED  STRESS  FROM  OBSERVED  VALUES 
FROM  SPL  TO  SI  371.  PSI  •••••  FROM  SI  TO  SMAX  38.  PSI 
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COMPARISON  OF  COMPUTEO  STRESS  VALUES 
FOR  SPECIMEN  NO.  ESROll 


WITH  EXPERIMENTALLY  OBSERVED  VALUES 
DATE  OF  TEST  OEC.  1974 


STRESS-ENGINEERING  STRAIN  STRESS-PLASTIC  STRAIN  COMPONENT 


OBSERVEO 

OBSERVED 

ADJUSTEO 

COMPUTED 

STRESS 

PLASTIC 

COMPUTED 

STRESS 

STRESS 

STRAIN 

STRAIN 

STRESS 

DEVIATION 

STRAIN 

STRESS 

DEVIATION 

KSI 

PCT 

PCT 

KSI 

KSI 

PCT 

KSI 

KSI 

• 00 

• 0C01 

.0C05 

.14 

-.140 

.CCGC 

• 00 

.000 

10.30 

.0342 

• 0346 

10.23 

.067 

• 0000 

.00 

.000 

2C.60 

.0699 

.0703 

20.80 

-.201 

.0000 

.00 

•COO 

30.90 

.1035 

.1039 

30.75 

.153 

.0000 

• 00 

.000 

41.20 

.1397 

.1401 

41.46 

-.262 

.0000 

•oc 

• 000 

51.50 

.1727 

.1731 

51.23 

.27C 

.COCO 

• 00 

•000 

61.30 

.2072 

.2076 

61.44 

.358 

.0000 

.00 

.000 

72. 1C 

.2433 

.2437 

72.13 

-.028 

.0000 

• 00 

.coo 

82.40 

.2772 

.2776 

82.16 

.237 

• 0000 

.00 

• oco 

92.70 

• 3128 

.3132 

92.70 

-.C01 

• ccoo 

• 00 

•ccc 

103.00 

• 3486 

.3490 

103.30 

-.298 

.0000 

•cc 

• OCO 

113.20 

• 3835 

.3839 

113.63 

-.428 

• cooc 

.00 

.ooc 

123.60 

• 4163 

.4172 

123.49 

.115 

.0000 

• 00 

.000 

133.8C 

.4512 

.4516 

133.67 

.132 

• ccoo 

.CO 

.occ 

144.10 

• 4865 

.4869 

144.12 

-.017 

.0000 

.00 

• 000 

154.40 

.5211 

.5215 

154.36 

.041 

.0000 

.00 

•ooc 

164.70 

.5562 

.5566 

164.66 

.044 

.0000 

.00 

• 000 

175.00 

.6021 

• 6025 

175.20 

-.199 

• 0113 

175.27 

-.266 

181.20 

.6701 

.6705 

180.62 

.582 

.0583 

180.56 

• 640 

183. 3C 

.7590 

.7594 

183.31 

- .006 

• 1401 

183.45 

-.148 

185.30 

.8863 

• 8867 

185.31 

-.010 

.2607 

185.52 

-.221 

187.40 

1.0471 

1.0475 

186.86 

.544 

.4144 

187.42 

-.024 

189.50 

1.2790 

1.2794 

188.71 

.788 

.6392 

189.44 

.055 

19C.3C 

1.4084 

1.4088 

189.61 

• 687 

.7659 

190.29 

.005 

190.70 

1.4890 

1.4894 

190.14 

• 565 

• 8451 

190.76 

-.059 

191.10 

1.5520 

1.5524 

190.53 

.575 

• 9068 

191.09 

.009 

TANGENT  MODULUS  VALUES  AND  POISSON’S  RATIO 
SPECIMEN  NUMBER  ESROll  OATE  OF  TEST  DEC.  1974 


COMPUTED 

STRESS 

KSI 


• 14 

10.23 
20.80 
30.75 
41.46 

51.23 
61.44 
72.13 
82.16 
92.70 

103.30 
113.63 
123.49 
133.67 
144.12 
154.36 
164.66 
175.20 
180.62 

183.31 

185.31 
186.86 
188.71 
189.61 
190.14 
190.53 


ADJUSTED 
STRAIN 
PC  T 


• 000 

• 035 
.070 

• 104 
.140 
.173 

• 208 
.244 
.278 
.313 
.349 
• 384 
.417 
.452 
.487 
.521 
.557 
.602 
.670 
.759 
.887 

1.047 

1.279 

1.409 

1.489 

1.552 


TANGENT 

MOOULUS 

PSI 


.2960*08 

•2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2960*08 

.2870*08 

.1464*08 

.4458*07 

.2127*07 

.1195*07 

•8818*06 

.7292*06 

.6653*06 

.6311*06 

.6067*06 


TANGENT 

MODULUS 

RATIO 


1.0000 

l.OCCC 

1.0000 

l.OCCC 

1.0000 

l.OCOC 

1.0000 

l.CCCO 

1.0000 

l.OCCC 

1.0000 

l.CCCO 

1.0000 

l.OCCC 

1.0000 

l.OCOO 

.9695 

• 4 S46 
.1506 

• C 71  9 
.0404 
.0298 
.0246 

• C225 
.0213 
.0205 


COMPUTED 

POISSON’S 

RATIO 


• 280C 

• 2 80C 

• 28C0 

• 2800 
.2800 
• 2800 
.2900 

• 280C 
.2800 

• 28C0 
.2800 
.2800 
.2800 
.2800 
.2800 
.2800 
.2801 

• 2839 
.2998 
.3206 

• 3447 
.3674 
.3904 

• 4 COO 

• 4C  5 1 

• 4 C88 


20 


STRAIN-HARDENING  COEFFICIENTS  AND  EXPONENTS  FROM  A TWO  EQUATION  ANALYSIS  OF  THE 
STRESS-PLASTIC  STRAIN  COMPONENT  DATA  USING  A VARIABLE  INTERMEDIATE  STRESS  VALUE. 

SPECIMEN  NUMBER  ESROX1  DATE  OF  TEST  DEC.  1974 


EQUATION  1 PARAMETERS  ARE  FOR  VALUES  OF  STRESS  FROM  175.00  KSI  TO  AN  INTERMEDIATE  POINT  IN  THE  R-  0 REGION 
EQUATION  2 PARAMETERS  ARE  FOR  VALUES  OF  STRESS  FROM  THE  INTERMEDIATE  POINT  TC  A MAXIMUM  STRESS  OF  191.10  KSI 


EQUATION  1 PARAMETERS  INTERMEDIATE  POINT  EQUATION  2 PARAMETERS 


OBSERVED  PLASTIC  SECANT 


COEFF. 

EXPONENT 

DEVIATION 

STRESS 

STRAIN 

MOOULUS 

COEFF. 

EXPONENT 

DEVIATION 

A1 

Q 

KSI 

KSI 

PCT 

RATIO 

A2 

R 

KSI 

, 2077*06 

•01874 

• 223 

183.30 

• 4144 

.815 

•2146406 

.02468 

• 014 

2066406 

•01810 

• 186 

185.30 

• 6392 

.706 

.2147406 

.02475 

.017 

,2071406 

•01840 

.158 

187.40 

• 0000 

• 6C4 

.2135406 

.02362 

.021 

2080406 

•01897 

• 165 

189.50 

• 0000 

.500 

.2145406 

.024  57 

.019 
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